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Abstract 

We consider the plane 3 body problem with 2 of the masses small. 
Periodic solutions with near collisions of small bodies were named by 
Poincare second species periodic solutions. Such solutions shadow chains 
of collision orbits of 2 uncoupled Kepler problems. Poincare only sketched 
the proof of the existence of second species solutions. Rigorous proofs 
appeared much later and only for the restricted 3 body problem. We 
develop a variational approach to the existence of second species periodic 
solutions for the nonrestricted 3 body problem. As an application, we give 
a rigorous proof of the existence of a class of second species solutions. 

1 Introduction 



Consider the plane 3-body problem with masses mi, 7712,7713. Suppose that 7713 
is much larger than r7ii, 7712, i.e. fi = (riii + m2)/m3 is a small parameter. Then 

"I1/7713 — /itti, 7712/7713 — ^a2, ai + a2 = 1. 

Let (71, ^2 G be positions of mi, m2 with respect to m3 (Poincarc's heliocentric 
coordinates) and pi,P2 G their scaled momenta. The motion of mi, m2 with 
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respect to is described by a Hamiltonian system (-ff^) with Hamiltonian 

H,i,,p) ^ Hoi,,p) + , (^^^ - , (1.1) 

V 2 \q1-q2\J 



where 



2ai 2a2 \qi\ \q2\ 

The Hamiltonian Hf^ is a small perturbation of the Hamiltonian Hq = Hi + 
H2 describing 2 uncoupled Kepler problems. The configuration space of system 
(Hq) is — U X U, where [/ = \ {0}. The configuration space of the 
perturbed system (-ff^) is \ A, where 

A = {q = (gi,g2) eU^ -.qi^ 52} 

represents collisions of rrii, 7712. 

System (-ff^) has energy integral and the integral of angular momenturrQ 

G{q,p) = Gi + G2 = iqi ■ pi + iq2 ■ P2 

corresponding to the rotational symmetry e'^ : — > M^. 

System (Hq) has additional first integrals Hi,H2 and Gi, G2 - energies and 
angular momenta of mi, m2. In the domain 

P = {{q,p) eU^xR^: Hi,H2 < 0, Gi,G2 ^ 0}, 

orbits of mi, m2 are Kepler ellipses, and solutions of system (Hq) are quasiperi- 
odic with 2 frequencies. 

Let R C P he the regular domain - the set of points in P such that the 
corresponding Kepler ellipses do not cross. For small fi > solutions of system 
(Hfj^) in R are 0(/i)-approximated by solutions of system (Hq) on finite time 
intervals independent of fi. By the classical perturbation theory, away from 
resonances the same is true on longer time intervals. Moreover, as proved by 
Arnold [3^, for small fi > 0, system (-ff^) has a large measure of invariant 
2-dimensional KAM tori on which solutions are quasiperiodic, and thus well 
approximated (modulo rotation) by solutions of system (Hq) on infinite time 
intervals. 

In the singular domain S G P, where the corresponding Kepler ellipses cross, 
the classical perturbation theory does not work. Indeed, for almost any initial 
condition in S, solution of system (Hq) is quasiperiodic with incommensurable 
frequencies, and so eventually mi,m2 simultaneously approach an intersection 
point of Kepler ellipses. Then the perturbation in (jl.ip becomes large, and so 
it can not be ignored even in the first approximation in fj,. 

For small ^ > solutions of the 3 body problem (7?^) in S can be described 
as follows. The bodies toi, m2 move along nearly Kepler ellipses and after many 
revolutions they almost collide. Then they start moving along a new pair of 



^We identify with C, so multiplication by i is rotation by 7r/2. 
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nearly Kepler orbits. If the new energies Hi,H2 are both negative, so the 
new Kepler orbits are ellipses, then mi, 1712 will again nearly collide after many 
revolutions, and the process repeats itself. Thus almost collision solutions of 
the 3-body problem (i?^) shadow chains of collision orbits of system (-ffo)- 

Almost collision periodic solutions of system (iJ^) were first studied by 
Poincare in New Methods of Celestial Mechanics. Poincare named them sec- 
ond species periodic solutions. However, he did not provide a rigorous existence 
proof. Rigorous proofs appeared much later (see e.g. [inilinilZ]) and only for the 
restricted 3 body problem, circular and elliptic. In [TUJ |S] also chaotic second 
species solutions of the circular and elliptic restricted problem were studied. 

The goal of this paper is to develop a variational approach to almost collision 
periodic orbits of the nonrestricted 3 body problem. As an application, we will 
give a rigorous proof of the existence of a class of almost collision periodic orbits. 
Chaotic almost collision orbits will be studied in another paper. 

Remark 1.1. It is possible to fix the value of angular momentum G and reduce 
rotational symmetry. Then we obtain a Hamiltonian system with 3 degrees 
of freedom. However, since reduction of the rotational symmetry considerably 
complicates the Hamiltonian, it is simpler to work with the original Hamiltonian 
system (H^) with 4 degrees of freedom. 

Remark 1.2. We consider only near collisions of small masses mi and m2 
and exclude near collisions 0/ mi, 7712 with m^. In particular, triple collisions 
are excluded. It is well known that double collisions can be regularized, but 
the Levi-Civita regularization becomes singular as pL Q. Understanding this 
singularity is the base for our methods. Levi-Civita regularization was previously 
used to study second species solutions for the restricted 3 body problem, see e.g. 

m\MM- 

Remark 1.3. Main results of this paper hold for more general Hamiltonians 
with singularity, for example 

Hf.{q,p) ^ 7r-r^ + Tr~r^ + 9[(i.iJ-)- \ r, ai,a2,/>o, 

2ai[q) 2a2(q) m - q2\ 

where all functions are smooth without singularity at qi — q2. 



2 Main results 

A solution of the Hamiltonian system (iJ^t) is determined by its projection to 

the configuration space f/^ \ A which will be called a trajectory. Let be the 
Lagrangian corresponding to iJ^. A T-periodic trajectory 7 : R — \ A is a 
critical point of the Hamilton action functional 

A^(T,7)= r L^{j{t),^{t))dt (2.1) 
-'0 
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on the space of T-periodic W^^^ curves in C/^ \ A. We write T explicitly in 
{T, 7) because later it will become a variable. 

Any trajectory of system (Hq) has the form 7 = (71,72), where 7^ is a 
trajectory of the Kepler problem. For ^ — the action functional A — Aq 
splits: 

A(T,7) =aiB(r,7i) + a2B(r,72), 7 = (71,72), (2.2) 

where 

is the action functional of the Kepler problem on the space of T-periodic W^^^ 
curves a -.R^ U. 

We have to consider also trajectories of system (Hq) with collisions. A T- 
periodic curve 7 = (71,72) : K — is called a periodic n- collision chain if 
there exist time moments 

t = (ti,...,t„), h < ... <tn <tn+l=h+T, (2.4) 

such that: 

• 7 has collisions at i = tj , so that 

li^j) = (xj-.xj) e A, 71 (*j) = l2{tj) = Xj. 



• 7l[t ,t +1] is a trajectory of system (Hq) which will be called a collision 
orbit. 

• Momentum p(t) = {pi(t),p2{t)) — (01171 (t), 0272 (t)) changes at collisions 
so that the total momentum y = pi + P2 is continuous: 

y{t,+0) = y{tj-0)=yj. (2.5) 

Equivalently, the jump of momentum p{tj + 0) — p{tj — 0) is orthogonal 
to A at 7(tj). 

• The total energy Hq does not change at collision: 

Ha{l{t,),p{t, + 0)) = Hoh{tj),p{tj - 0)). (2.6) 

By (|2.5p . the total angular momentum G = G1 + G2 is preserved at collisions: 
G'(7(^j),P(^j ±0)) =ixj -yj. 

Hence G is constant along 7. By (|2.6p . the total energy Hq = Hi + H2 = E \s 
also constant along 7, but not the energies Hi,H2 of mi,m2, or their angular 
momenta Gi , G2 . 

A collision chain 7 is a broken trajectory of system (Hq) ~ a concatenation 
of collision orbits with reflections from A. However, unlike for ordinary billiard 
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systems, A has codimension 2 in the configuration space U'^, so the change of the 
normal component of the momentum at coUision is not uniquely determined. 
Thus there is no direct interpretation of collision chains as trajectories of a 
dynamical system. Such an interpretation is given later on. 

Collision chains are limits of trajectories of system (-ff^t) which approach 
collisions as fi 0. Indeed, we have: 

Proposition 2.1. Let 7^ be a Tf^-periodic trajectory of system {H^) which 
uniformly converges, as fi 0, to a T-periodic curve 7. //7([0, T]) H A is a 
finite set, then j is a periodic collision chain. 

We say that 7^ is an almost collision orbit shadowing the collision chain 7. 
A similar statement holds for nonpcriodic collision chains. 

Intuitively, Proposition 12.11 is almost evident: a near collision of mi, 7712 
lasts a short time during which the influence of the non-colliding body is 
negligible. Then mi,m2 form a 2 body problem, so their total momentum 
y = Pi -\- P2 and total energy Hq = Hi + H2 are almost preserved. This yields 
(|2.5I) - (|2.6I) . This can be made into a rigorous proof, see e.g. T. A better way 
to prove Proposition 12 . 1 1 is by using the Levi-Civita regularization, see section 



Collision chains can be characterized as extremals of Hamilton's action func- 
tional (12.21) . Let be the set of a; = (t,T, 7), where t = {ti, . . . ,t„) satisfies 
(|2.4|) and 7 : R — > [/^ is a T-periodic W^^^ curve such that j{tj) = {xj,Xj) G A. 
The collision times tj and collision points Xj are not fixed. Then fij can be 
identified with an open set in a Hilbert space (see (|2.8p and collision times tj 
and collision points smooth functions on Vli . 

Remark 2.1. If ^{t) ^ A for t 7^ tj, then time moments tj are determined 
by the curve 7, i.e. the projection W^''^{M./TZ,U'^), (t,r, 7) — > 7, is 

injective. Then lo — (t,T, 7) is determined by 7. But tj are not continuous 
functions ofj, so we have to include the variables t in the definition offl^. 

Remark 2.2. In the one- dimensional calculus of variations Hilbert spaces are 
unnecessary: at least locally function spaces can be replaced by finite dimen- 
sional subspaces of broken extremals. We will use this approximation later on. 
However, in this section we use conventional W^''^ setting. 

The action functional A{uj) = A(T, 7) is a smooth function on il^. 

Proposition 2.2. j is a T-periodic n-collision chain iff lu = (t,T, 7) is a 
critical point of the Hamilton action A on il^. 

Proof. If w = (t, T, 7) is a critical point of on fl^, then each segment 7|[ti_i,tj] 
is a collision orbit of system (Hq). Then by the first variation formula [J, 



El 



□ 



n 
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where A/ij and Aj/j are jumps of energy and total momentum: 

Ay,^y+-yJ, yf^y{t,±0), Ah,^h+-hJ, hf ^ Ho{l{t,),p{t,±Q)). 

Since the differentials dtj , dxj are independent, critical points of A satisfy Ahj = 
and Ayj = 0. This implies p.5p - (|2.6p . Converse is also evident. □ 

For collision chains with fixed energy Hq — E we use Maupertuis's varia- 
tional principle [2]. Let fin — Ut>o5^^- The period T > is now a smooth 
function on Hamilton's action is replaced by the Maupertuis action func- 
tional 

^^H = A^(r,7) = A(T,7) + i?T, ^ = (t,r,7). (2.7) 
If 7 is a collision chain with energy E, then 

A^(r,7)= fp-dq 

is the classical Maupertuis action. The Maupertuis principle for collision chains 
is as follows: 

Proposition 2.3. 7 is a T -periodic n-collision chain of energy E iff ui = 
(t,T, 7) is an extremal of the Maupertuis functional on f2„. 

Due to time and rotation invariance, critical points of the action functional 
are degenerate. To eliminate time degeneracy, we identify collision chains which 
differ by time translation 

7(i)^7(i-T), tj^tj+T, T e K. 

This defines a group action of M on . The coordinates on the quotient space 
/K can be defined as follows. Let 

Sj ^tj+i-tj > 0, ajir) =-/j{tj+TSj), 0<T<1. 

Denote 

s = (si, . .. ,s„) e M!;:, cr = (cti,. . . ,a„) e W^, 

where 

Wn = {a^iai,.. . ,a„) : a, G W^'\[0, 1], U^), a,il) = fT,+i(0) e A}. 

The map map ri„ — ^ K" x Wn, (t,r, 7) (s, cr), makes it possible to identify 
i7„ with M" X Wn- We can represent aj by 

a, e <''([0, 1], C/2), a,(r) = ^.(t) - (1 - r)x, - tx,+i. 

Then a is determined by (x, a), where x = . . . , x„) € C/" and a = (cti, . . . , an)- 
Hence Wn = C/" x Wq '^([0, 1], t/^"). Finally 

hn=RlxWn = M.lxU'' xWo^^{[0,l],U^'^), = f2„ X R. (2.8) 
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The action functional gives a smooth function A{s, a) on f2„, invariant under 
rotations: ^(s, e^^a) = A{s, a). We deal with the rotation degeneracy later on. 

Often it is convenient to use parametrization independent Jacobi's form of 
the Maupertuis action functional - the length of 7 in the Jacobi's metric dss- 

J^il) = / dsE, dsE = max{p • dq : Ho{q,p) = E} 

J-y P 

= V2(£^ + ai|gi|-i +a2|52|-i)(Q;i|(i5i|2 + a2|rfg2|2). 

Then J^{j) < A^{T,^) and if 7 is parametrized so that Hq = E, then 
A^(T, 7) = J^(7). Thus up to parametrization, extremals of and are 
the same. For a collision chain 7 corresponding to (s,0-) G f2„, 

n 

is a function on Wn- We obtain 

Proposition 2.4. If ^ is a periodic n-collision chain of energy E then the cor- 
responding u e Wn is an extremal of the Jacobi action . If a is an extremal 
of , then if each Gj is reparametrized so that Hq = E, the corresponding 7 

is a periodic n-collision chain. 

We will not use Jacobi's variational principle since is not a smooth 
function. A discrete version of Jacobi's action is smooth and we will use it later 
on. 

Similar variational principles hold for collision chains periodic in a rotating 
coordinate frame: 7(4 + T) = e'*7(f) for some quasiperiod T and phase We 
call 7 periodic modulo rotation. If $ ^ 27rQ, then 7 is quasiperiodic in a fixed 
coordinate frame. 

The corresponding function space is defined as follows. Let n„ be the set of 
all (t,T,7, $), where t,T arc as before, $ e R and the W^^^ curve 7 : R 17^ 
satisfies 7(4 + T) = e**7(t) and ^{tj) e A. Then Qn can be identified with an 
open set in a Hilbert space and tj,Xj,T, $ are smooth functions on n„. In fact 
f2„ ^ f2„xIR. Indeed, j{t) = e-***/'^7(t) is a T-periodic curve, so {t,T,j) e n„. 

Define the Maupertuis-Routh action functional on n„ by 

>l^°(t,T,7,$) = A(T,7) + £;T-G$. (2.9) 
This is a smooth function on n„ and we have: 

Proposition 2.5. 7 is a periodic modulo rotation collision chain with energy 
E and angular momentum G iff (t, T, 7, $) is a critical point of the functional 
A^'^ on Vln. 

Remark 2.3. It seems natural to take $ € T = M/27rZ since $ + 27r gives the 
same collision chain. But then the functional A^'^ will be multivalued: defined 
modulo 27rG. 
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Remark 2.4. The name of the functional AP'^ is motivated as follows. One 
can perform Routh 's reduction of the rotational symmetry for fixed G replac- 
ing the configuration space U"^ by — U"^ /T = x T and the Lagrangian 
by the so called Routh function. Then the functional A^'~^ becomes the Mau- 
pertuis functional for the reduced Routh system. Probably this observation is 
due to Birkhoff /5/. However, Routh's reduction makes the Lagrangian more 
complicated, so we do not use it in this paper. 

There are several other possible variational principles for collision chains (for 
example, we may fix the phase $), but in the present paper we will use only the 
ones given above. 

SufScient condition for the existence of a periodic orbit of system (ff^ ) shad- 
owing a given collision chain 7 requires that the chain is nontrivial in the fol- 
lowing sense. Let u{t) — 72(t) — 7i(i) and let 

v(t) = au{t) — aip2{t) — a2Pi{t), a = a\a-2., (2-10) 

be the scaled relative velocity of mi, TO2. Let = v{tj ± 0) be relative collision 
velocities. Since 

/.± = |y,f/2+|«±|V2a-|x, (2.11) 

equations (|2.5p - p.6|) imply that = \u'^\'. relative speed is preserved at 
collision. We impose two essentially equivalent conditions: 

Direction change condition. Relative collision velocity changes direction at 
collision: Aw^ — — v~ ^ 0, j — I, . . . , n. In particular, ^ 0. 

No early collisions condition. 7(t) ^ A for t ^tj. 

If the direction change condition is not satisfied at some tj, then 7(tj — 0) = 
7(tj -|-0), and so 7|[t ._^ j .^^j is a smooth trajectory of system {Hq). Deleting the 
collision time moment tj we obtain a (n — l)-collision chain violating no early 
collisions condition. 

Conversely, if 7 is a 71-collision chain violating no early collisions condition, 
then adding an extra collision time moment, we obtain a (n -I- l)-collision chain 
violating the changing direction condition. From now on we add these two 
equivalent conditions to the definition of a collision chain. 

Remark 2.5. The changing direction condition implies that almost collision 
orbits 7^ shadowing the collision chain 7 come 0{^)-close to collision. Often 
almost collision orbits discussed in Astronomy come close to collision, but not 
too close, for example 0{iJi'^) -close with v e (0, 1), see e.g. /J 51 [HI [7^ . Such 
orbits change direction at near collision, but this change is small as fi ^ 0. 
Then the corresponding collision chains do not satisfy the changing direction 
condition: Avj = v^'—vj = 0. Our methods do not work for such almost 
collision orbits. 
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The changing direction condition makes it possible to construct a shadowing 
orbit 7^ of system (-ff^), but it does not prevent 7^ from having regularizable 
double collisions of mi,m2- To exclude such collisions we need to impose an 
extra condition: 

No return condition, v'j' + vj ^ 0, j = 1, . . . ,n. 

But this condition is not as crucial as the changing direction condition, so 
we do not include it in the definition of a collision chain. 

To construct shadowing orbits we also need some nondegeneracy assump- 
tions. We say that a T-periodic 71-collision chain 7 with energy E is nondegen- 
erate if w = (t,T, 7) £ f2„ is a nondegenerate modulo symmetry critical point of 
the Maupertuis action on ri„. 

Due to time translation and rotation invariance, critical points of are all 
degenerate: the group action j{t) — >■ e'^^j{t — r) of M x T preserves A^ . We 
say that uj — (t,T, 7) e Qn is nondegenerate modulo symmetry if the nullity 
of the quadric form d'^A{uj) on T^^iln is 2 - the lowest possible. Equivalently, 
the manifold M C f^n obtained from uj by the action of the group M x T is a 
nondegenerate critical manifold. Nondegeneracy modulo symmetry is equivalent 
to nondegeneracy of the corresponding critical point on r2„/ (R x T) = f2„/T. 

As usual in the classical calculus of variations, the Hessian operator cor- 
responding to d'^A^{uj) is a sum of invertible and compact operators on the 
Hilbert space Tt^f2„, so nondegeneracy modulo symmetry implies that the Hes- 
sian has bounded inverse on T^^ln/T^^M . In fact, at least locally, A^ can be 
reduced to a finite dimensional discrete action functional (see section [5]), so all 
Hilbert spaces involved are essentially finite dimensional. 

Now two main results will be formulated. 

Theorem 2.1. Let ^ be a nondegenerate T-periodic collision chain with energy 
E. Then for small fi > there exists a T^j^-periodic orbit 7^ of system (H^) with 
energy E which 0{ii) shadows 7; 

T^^T + Oifi), 7mW=7W + 0(m), te[0,T]. 

If 7 satisfies the no return condition, then 7^ has no collisions and there exist 
< a < b independent of /i such that 

^j.a<d{j^{tj),A) <fib. (2.12) 

Due to time translation and rotation symmetry, 4 multipliers of 7^ (eigen- 
values of the linear symplectic Poincare map of R^) are equal to 1. Nontrivial 
multipliers are Ai, Aj^^, A2, A^^, where Ai(/^) is real and large of order | ln/x|, and 
A2(/i) has a limit as /i — > 0, where A2(0) ^ 1 is complex with |A2(0)| = 1 or real. 

Next we consider collision chains with fixed energy E and angular momentum 
G. Again we say that a periodic modulo rotation collision chain 7 is nondegener- 
ate if the corresponding (t, T, 7, $) G r2„ is a nondegenerate modulo symmetry 
critical point of the functional A^^ on f2„ . Thus it has only degeneracy coming 
from rotation and time translation invariance. 
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Theorem 2.2. Let j be a T -periodic modulo rotation nondegenerate collision 
chain with energy E and angular momentum G. Then for small fj, > there 
exists a periodic modulo rotation orbit 7^ of system {H^) with energy E and 
angular momentum G which 0{^)-shadows 7; 

7;.(i + r^) = e^*''7/.W, 7f.it) = lit) + 0{f^), te[0,T], 

where 

T^ = T + Oiti), $^ = $ + 0(//). 

An estimate (|2.12p holds also here. Even if 7 is periodic ($ S 27rQ), in 
general the shadowing orbit 7^ will be periodic only modulo rotation, and thus 
quasiperiodic in a fixed coordinate frame. 

To use Theorems I2.1H2.21 we need to find nondegenerate modulo symmetry 
collision chains. In general this is not easy. A simple application of Theorem 
12. 2[ based on a perturbative approach, is given in section [3] More complex 
applications will be given in a future publication. 

In section [5] a description of nondegenerate collision orbits is given. Using 
this description, in section [5] we reduce the action functionals to their discrete 
versions. Then in section [6] we formulate a local connection result - Theorem 
16.11 - and use it to prove Theorem 12. II The proof of Theorem 12.21 is similar. In 
section [7] we use Levi-Civita regularization to reduce Theorem 16.11 to Theorem 
I7.2l which is a generalization of the Shilnikov Lemma [H] to Hamiltonian systems 
with a normally hyperbolic critical manifold. 

Remark 2.6. In this paper we do not attempt to use global variational methods. 
The reason is that although one can use global methods to find critical points 
of the action functionals, in general it is hard to check that the critical points 
satisfy the changing direction condition. 

3 Restricted elliptic limit 

Suppose that one of the small masses mi, 7712 is much smaller than the other: 
ai <^ a2. In the formal limit ai — >■ we obtain the restricted elliptic 3 body 
problem for which many second species periodic solutions were obtained in [7]. 
These results do not immediately extend to the case of small ai > 0. However, 
we will show that they can be used to obtain many second species periodic 
solutions for the nonrestricted 3 body problem. 

Let us fix energy E and angular momentum G. For ai <^ a2, the Maupertuis- 
Routh action functional 

^^'=(t,T,7,$) = aiB(T,7i) + a2B(T,72) + ^r-G$, 7 -(71,72), (3.1) 

on Qn is a small perturbation of the Maupertuis-Routh action functional for 
the Kepler problem. Indeed, the functional Aq"-^ = A^'-^\ai=o does not depend 
on 7i: 

Ao'^it, T, 7, $) = B^^{T, 72, $) = B{T, 72) +TE- G$. (3.2) 
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The condition j{tj) G A imposes no restrictions on 72. Thus the functional B^^ 
is defined on the set 11 of (T, a, <&), where a : R — >■ i7 is a W^^^ curve such that 
a{t + T) = e**(7(t). We have a submersion tt : n„ ^ H, (t, T, 7, ^ (T, 72, $), 
and y^f ^5 = o 7T. 

The functional S^*^ is very degenerate, because all orbits of the Kepler 
problem with energy E <0 are periodic with the same period r = 27r(— 2i?)~^/^. 
Suppose E, G arc such that there exists an elliptic orbit F of Kepler's problem 
with energy E and angular momentum G. For definiteness let G > 0. Then 
< {—2E)G < 1. The major semiaxis and eccentricity of F are 

a={-2E)-\ e = Vl + 2EG. 

The Maupertuis action is 

Je{T) = j ydx = 2'K{-2E)-^''^. 

The counterclockwise elliptic orbit F : R — ^ /7 is defined uniquely modulo rota- 
tion and time translation. 

Proposition 3.1. Let E <(), G > and {—2E)G < 1. Then all critical points 
uj = {T, a, of the functional B^"-^ on H belong to one of the nondegenerate 
critical manifolds Mm C U, m £ N, obtained from {mr, F, 0) by rotation and 
time translation ofT. We have 

B^^\m^ = B{mT,r) + mrE = mJEiT) = 2nm{-2E)-^^'^ . 

Proof. Let (T, cr, <I>) G 11 be a critical point of B^"^ . Then cr is a solution of the 
Kepler problem with energy E and angular momentum G and hence cr is a time 
translation and rotation of F. Since F is a non-circular orbit, quasiperiodicity 
condition a{t + T) = e'*cr(f) implies that $ = mod 27rZ and T = mr for some 
m e N. 

Next we need to check that is a nondegenerate critical manifold of 
B^'^. Essentially this is the same statement, but now we need to consider the 
linearized Kepler problem. 

The second variation d^B^^{(jj) at w = (mr, F,0) is a bilinear form on 

the tangent space T^H which is the set oi rj = {6,£,,(j)), where 6,(j) £ M. and 
^ : M — )■ is a vector field such that 

^{t + mr) = ^{t) + t{t)9 + ir{t)(l). (3.3) 

The standard calculus of variations implies that if ry G T^H belongs to the 
kernel of d^B^'^{uj), then ^ is a solution of the variational equation for F which 
lie on the zero levels of the linear first integrals corresponding to the integrals of 
angular momentum and energy. The linear approximations at F to the integrals 
of energy and angular momentum are 

m ■ m - m ■ m = o, im ■ m - im ■ m = o. 
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Condition p.3p gives 

iT{t) ■ f{t)0 EE 0, iT{t) ■ f (t)0 EE 0. 

Since F is noncircular, 9 — (j> — and so ^(t + mr) = ^(t). It follows that ij = 
(0, 5, 0) is tangent to Mm, i-e. the variation ^(i) is obtained by time translation 
and rotation of r{t). □ 

The critical manifold Nm — 7i'~^(A/m) C (in of Aq'^ corresponding to Mm 
is (up to time translation and rotation) 

Nm = {(t,mT,cr,r,0) e Cln : cr(t + mr) = cr(tj) = r(ij)}. 

This is an infinite dimensional nondegenerate critical manifold of Aq'^ . For 
nonzero ai, by p.ip . 

A'^^In^ = ai(S(T, ct) + i;TOr - 27rTO(-2i;)-i/2) + 27rm(-2£;)-i/^ 

By a standard property of nondegenerate critical manifolds TT , any nonde- 
generate modulo symmetry critical point oj G Nm of A^'^\n,„ for small ai > 
gives a nondegenerate modulo symmetry critical point of A^'~^, and hence a 
nondegenerate modulo symmetry collision chain with energy E and angular 
momentum G. 

Up to an additive constant and a constant multiple, A^^In^ Hamilton's 
action B{mT, a) for the Kepler problem. It is defined on the set Ilr^m of (t, cr), 
where tr : K — >■ [/ is an mr-periodic curve such that <j{tj) = T{tj). Thus 
B[mT,a) — ,Br,m(t,c) is precisely the action functional whose critical points 
are collision chains of the elliptic restricted 3 body problem. This functional 
was studied in [7], and many its nondegenerate critical points were found for 
small eccentricity (almost circular F), i.e. {—2E)G close to 1. Also the changing 
direction and no early collisions condition was verified in [7j, and this carries 
out for small ai > 0. We obtain 

Theorem 3.1. Let < {~2E)\G\ < 1 be close to 1. Then for sufficiently small 
ai > there exist many collision chains 7 such that for sufficiently small fi > 0, 
7 is 0{fi) -shadowed by a second species periodic modulo rotation solution 7^ of 
the nonrestricted 3 body problem with given E, G. 

This result can be improved by using a more quantitative statement from 
[7]. The obtained second species solutions are periodic in a rotating coordinate 
frame and quasiperiodic in a fixed coordinate frame. Proper periodic orbits will 
be obtained in a future publication; for them reduction to the restricted elliptic 
problem is impossible. 

4 Collision action function 

Collision chains can be represented as critical points of a function of a finite 
number of variables ~ discrete action functional. This is needed for the proof 



12 



of Theorems UTTHUll and in subsequent publications. Since collision chains are 
concatenations of collision orbits, we need to describe collision orbits first. 

A collision orbit 7 — (71, 72) of system {Hq) is a pair of Kepler orbits joining 
the points x^,x^ e U. Thus description of collision orbits is reduced to the 
classical Lambert's problem [52] of joining the points x^,x^ by a Kepler orbit. 

First we join the points X-,x+ by a Kepler orbit F : [0,t] — J7 with fixed 
energy E < 0, or, equivalently, fixed major semi axis a — {~2E)~^. Due to 
scaling invariance of Kepler's problem without loss of generality set a = 1. 
Then a Kepler ellipse passing through x_ , x+ is determined by the second focus 
F such that 

|x_| + |a;_-F|-2, \x+\ + \x+ - F\ ^ 2. (4.1) 

The solution F = F{x-, x^) of these equations exists and smoothly depends on 
x± if the corresponding circles intersect transversely, i.e. {x^,x+) lie in the set 

X = {{x^,x+) e C/^ : ||a;+| - |a;_|| < |a;+ - a;_| < 4- - |a;+|}. 

For {x-tX+) € X there exist two solutions F of equations ()4.ip . and we take 
one of them, for definiteness the one on the left side of the segment x-Xj^. 

Let F(a;_, a;+) be the counter clock wise simple arc of the constructed Kepler 
ellipse joining the points a;_ and X-\-. Let 

f{x-.x+) ^ J^ydx = ^(2Nr' - l)'/^\dx\ 

be the Maupertuis action of F. This is a smooth rotation invariant function 
on X : 

f{e''^x-,e'''^x+) = f{x-,x+). 

Remark 4.1. By Lambert's Theorem J2^, f is a function of s± — |a;_| + |a:+|± 
|x_ — a;+| only. An explicit formula is 

fix_,x+) = W{s+)±W{s_), (4.2) 

where 

^(s) = 2 V(4 - s)s + 2 arctan 

Plus is taken if xj^- = e*^a;_ with 9 £ [tt, 27r) and minus if £ (0,7r]. One can 
check that f is smooth at 9 — tt. 

Due to scaling invariance of the Kepler problem, for arbitrary negative en- 
ergy i? < 0, the Maupertuis action of a simple counter clock wise arc F = 
F(£', x-, a;+) connecting the points (x_,a;+) e Xe — {—2E)~^X is 

f{E,x,,x+) = J^ydx^ jm^-^ + E)f'^\dx\ 
= {-2E)-'/''f{{-2E)x^,{-2E)x+). 
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Next we describe Kepler orbits F : [0,t] — ;> U connecting x-,x+ while 
making n = [T] full revolutions around 0. To define the number n G Z of 
revolutions, set 

n = [{0{t) - 9{0))/27t], T{t) = r{t)e'^^'K 
Proposition 4.1. For any (a;_,a;_|_) e Xe and any n E Z; 

• There exists a Kepler orbit T — r„(_E, x_ , a;+) : [0,t] U of energy E 
joining the points x_ , 0:4. and making [T] — n revolutions. 

• r smoothly depends on _E, a;_,x+. 

• The Maupertuis action ofT is 

Jn{E,x^,x+) ^ Jy-dx = 
= {-2E)-^/'^{2Tr\n\ + {sgnn)f{{-2E)x^, {-2E)x+)). (4.3) 



d 

T = Tn{E,X-,X+) = —Jn{E,X_,X+). (4.4) 

Jn{E,x-,x+) > 0, {x-,x+)£Xe. (4.5) 



dE^ 



The orbits T„{E, x-, x+) are nondegenerate (have non-conjugate end points) 
and any nondegenerate connecting orbit with £' < is obtained in this way. 

For the classical Lambert's problem [22], when F is a simple elliptic arc, 
n = Oorn = — 1 depending on if F is a counterclockwise or clockwise. We set 
sgnO = 1. 

The first term in (|4.3|) is the Maupertuis action for n complete revolutions 
around the Kepler ellipse, and the second is the action of a simple elliptic arc. 
Equation (j4.4p follows from the first variation formula; it is essentially Kepler's 
time equation. So only inequality (j4.5p is non-evident. It is enough to check it 
for the classical Lambert's problem with n = 0, 1. Then (|4.5I) can be deduced 
from the explicit formula (14. 6p . although the computation is not trivial. An 
equivalent statement was proved in [19J. □ 

Next we consider Lambert's problem for fixed time r > 0. This problem 
involves solving the transcendental Kepler's equation so there is no explicit 
formula for the solution. Let D„ C R+ x U"^ be the open set which is the image 
of the diffcomorphism 

{E,x-,x+) ^ {Tn{E,x^,x+),x^,x+), E<0, {x^,x+)eXE. 

Proposition 4.2. For any n eZ and any {t,X-,x+) G Z3„; 

• There exists a nondegenerate Kepler orbit F = F„(t, a::_ , a;+) : [0,r] — )• U 
with [F] = n full revolutions joining the points x^ and X-j-. 
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• r smoothly depends on r, 

• Hamilton's action 

B{t,T) ^ F^{t,X-,x+) 
is a smooth function on Dn and 

— ^F„(t,x_,x+) < 0, (T,x_,a;+) e £>„. (4.6) 

• All nondegenerate connecting orbits with E < are T„ {t, 2;+) for some 
ri G Z and {t,x^,x+) € Dn- 

Proof. We need to find the energy E < such that the connecting orbit F = 
r„(i?, a;_, x+) in Proposition 14. II has given time r = t„{E,x^^x^). Then 

{F„{t,x^,x+) + tE) = J„{E,x-,x+). 

Hence J„ and — are Legendre transforms of each other: 
Jn{E, X- , = max(_F'„(r, a;_, + t£'), 

T 

i^„(T, a;_,a;+) = min( J„(i?, a;_, a;+) - tE). 

E 

Since J„(-E, x+) is convex with respect to E^ its Legendre transform —Fn{T, a;_|_) 
is convex in r and smooth. □ 

The initial and final total momenta of 7 are given by the first variation 
formula 

d d 
y+ = -Q^Fn{T,x-,x+), = F„(T,x_,a;+). (4.7) 

Now it is easy to describe nondegenerate collision orbits 7 = (71,72) of 
system (Hq). Denote by A: = [7] = (fci, ^2) S kj = [jj], the rotation vector 
of 7. We obtain 

Proposition 4.3. For any k = (fci,fc2) G Z'^ and any (r, a;_,a;+) G = 
Dk^ n Dk^ : 

• There exists a nondegenerate collision orbit j : [0,t] — C/^, 7 = 7(fc, r, a;_, a;+), 
wif/i collision points 7(0) = (a;_,a:;_), 7(r) = (x^,x+). 

• 7 smoothly depends on (r, a;_,a;+) G Vfc. 

• Hamilton's action of 7 is 

Sk{T,X-,x+) = A(r,7) = aiFk^{T,X-,x+) + a2^fe2 ("T: , a;+). (4.8) 



-^Sk{,T,x_,x+) <0, (r, a;_,x+) G Vfe. (4.9) 
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By the first variation formula [2], 

dSk{T,X-,x+) ~ y+ ■ dx+ — ■ dx^ — E dr, 

where E is the total energy of the collision orbit 7, and y± — y{t±) are total 
momenta at collisions. Thus 

y+ ^ -Q^[T,X-,x+), y_ (T,a;_,a;+), E = — —{t,x^,x+). 

(4.10) 

Remark 4.2. We will not need this in the present paper, but for almost all 
{t, x^ , x+) G Vk the collision action satisfies the twist condition 

det-^^{r,x^,x+)^0. (4.11) 
ox^ox+ 

Thus Sk is the generating function of a symplectic collision map (r, x_,?/_) — > 
(t, a;+,y+). This will be important for the study of chaotic collision chains. 

Remark 4.3. It can happen that the collision orbit ^ has early collisions: "/{t) G 
A for some t £ (0,t). To avoid this, we may need to delete from Vk a zero 
measure set, see 

Let us now fix energy E < {) and look for collision orbits of system [Hq) 
with energy E. The map 

{t,X-,x+) (-^^(t,x_,x+),x_,x+) 

is a diffeomorphism of Vk onto an open set Wk C M x C/^. Let be the 
Legendre transform of —Sk with respect to t: 

Lk{x-,x+) = max(S'fc(r,a:_,x+) -^tE) = {Sk{T,x-,x+) + tE)\^^^^(e,x^,x+), 

where r is obtained by solving the last equation (I4.10p . Then is a smooth 
function on 

= {(x_,x+) : (E,x^,x+) GWk}. 

We obtain 

Proposition 4.4. Let E < 0. For any {x^,x^) G there exists a unique 
collision orbit 7 : [0,r] U"^ of energy E such that 7(0) = (x_,a;_), 7(t) = 
(x+,a;+). Its Maupertuis action 

Lk{x-,x+)^ I p-dq ^ I dsE 
is a smooth function on . The total momenta at collision are 

^%x-,x+), y.^-^{x-,x+). (4.12) 



y+ 



dx+ ' ' dx 



16 



Ill terms of actions functions (|4.3p for the Kepler problem, 

Lf(a;„,a;+)= min (ai Jfe^ (£'1, a;_, a;+) + aa^/fe (-£^2, 

Remark 4.4. _Dwe to homogeneity of the Kepler problem^ 
where Lk corresponds to energy E = —1/2. 

Remark 4.5. The action functions Sk and can not be expressed in elemen- 
tary functions. However they admit simple asymptotic representation for large 
k. This will be done in a subsequent publication. 

In the next section we use the action functions Sk and to represent 
collision chains as critical points of discrete action functionals. 

5 Discrete variational principles 

For a given sequence k = (fc^, . . . , fc") e Z^", y e Z^, define a discrete Hamil- 
ton's action by 

n 

^k(s,x) ^^Ski{Sj,Xj,Xj+i), 

where s = (si, . . . , s„) G M" , x = (xi, . . . , Xn) & and S^j is the action 
function on V^i defined by (14.81) . The domain of is 

Vk = {(s,x) e W\_ X C/" : (sj,Xj,Xj+i) e Vfc, , a;„+i = a;i}. 

Any (s,x) e Vk defines (t,T, 7) g r2„ as follows. Take t = (ti,...,i„) so 
that Sj —tj+i — tj and set 

-i(t) ^ -f{y ,sj,xj,xj+i){t~tj), tj <t<tj+i, 

where -f{y ,Sj,Xj,Xj+i) : [0,Sj] — is the collision orbit in Proposition 14.31 
Then 7 = 7k(s,x) is a broken trajectory of system {Hq) with period T = 
X]J=i Hamilton's action ^k(s,x) = ^(T, 7). Of course (t,7) is defined 

modulo time translation, so we identify curves which differ by time translation. 
Thus we defined an embedding i : Vk — — f^n/IR and Ak ~ Ao i. 
For collision chains with fixed period T, we restrict to 

n 

yk^ = {(s,x)eyk:5].s,=T}. 

Proposition 5.1. Any critical point (s,x) 0/ Ak on defines a T-periodic 
collision chain 7 — 7k(s,x). 
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Indeed, critical points of satisfy 

-{Sj,Xj,Xj+i) H [Sj-i,Xj-i,x,,) = 0, (5.1) 



dxj ^"^'-^'"-^-^^ ' dXj 
dS 

-^{sj,Xj,Xj+i) = -E, (5.2) 

where —E is the Lagrange multiplier. By (|4.10l) . E is the energy of the corre- 
sponding collision chain 7. The total momentum at collision is 

Thus 7 satisfies □ 

Proposition 15.11 follows also from Proposition 12.21 Indeed, the functional 
is the restriction of Hamilton's action A to the set i(Vk) C i7„ of broken 

extremals. This set is obtained by equating to zero the differential of A for fixed 

t,T,x. 

In Proposition l5. li the period T is fixed. For collision chains with fixed energy 
i? < we consider a discrete Maupertuis action functional on Vk: 

N 

A^(s,x) = Ak(s,x) + i;r, T = Y.^i- 



Now T is a function on Vk, and A^(s,x) is the Maupertuis action (|2.7p of the 



broken trajectory 7 = 7k (s, x). We obtain 

Proposition 5.2. To any critical point (s,x) of on Vk there corresponds a 
periodic collision chain 7 with energy E. All nondegenerate collision chains with 
energy E are obtained in this way from nondegenerate modulo rotation critical 
points of some A^ . 

Remark 5.1. Hamilton's action is invariant under rotations: Ak(s,e'^x) = 
^k(s,x). Thus every critical point of the functional A^ is degenerate. To 
obtain nondegenerate critical points we should consider the quotient functional 
A^ on the quotient space 

Vk = Vk/T c M" X c7", [7" = [/"/T = wix r"-\ 

Let us now fix energy E < and angular momentum G and consider periodic 
modulo rotation collision chains 7 with given E, G. We obtain the discrete 
Maupertuis-Routh action functional 



^f^(s,x,$) = ^5fc.(sj,x„a;,+i)+^;r-G$, x^+i = e'*xi, ^ = 



AT 



The independent variables are s = (si, . . . , s„), x = (xi, . . . , x„) and <&, so the 
domain of A^'^ is 

Vk = {(s,x,$) e M!^ X [/" X R : (s,, a;,, x^+i) G Vy , x„+i = e'^xi}. 
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Proposition 5.3. To any critical point (s, x, $) of there corresponds a pe- 
riodic modulo rotation collision chain 7 = 7k (s, x, $) with energy E and angular 
momentum G. Any nondegenerate periodic modulo rotation collision chain with 
energy E and angular momentum G is obtained from a nondegenerate modulo 
rotation critical point of some A^'-^ . 

To construct orbits of system {H^) shadowing the colhsion chain 7 corre- 
sponding to a critical point (s,x), we need to verify the changing direction 
condition. For k = (fci, k^) € denote 

i?fc(T,a;_,x+) = Ffci(r, a;_,x+) - Ffe^ (r, a;_, 

By (|4.7p the relative collision velocities (|2.10p of a collision orbit 7 = 7(A:, r, a;_ , a;+) : 
[0, r] [/^ are given by 

■"(0) = --^{t,x-,x+), uij) = -^—{t,x-,x+). 

Thus the changing direction condition for the collision chain corresponding to 
(s, x) can be expressed as follows: 

-^{s-j,Xj,Xj+i) + —^—{sj-i,Xj-i,Xj) / 0. (5.3) 

We have 

Ak = aiBki + a2-Bk2, (5.4) 
where k = (ki, ka) with k^ = (fcj, . . . , kj) e Z" and 



Buj (s, x) = ^ (sj, x^, x^+i) (5.5) 

i=l 

is the discrete action functional for the Kepler problem. If (s, x) is a critical 
point of Ak with respect to x, then by (|5.ip . the changing direction condition 
(|5.3p is equivalent to 

^i3ki(s,x)=£0, j = l,...,n. (5.6) 



Next we reformulate the shadowing Theorems 12 . 11 - 12^ 

Theorem 5.1. Let (s,x) G Vk 6e a nondegenerate modulo rotation critical point 
of A^ satisfying the changing direction condition i5.6\} . Then for sufficiently 
small > the corresponding T-periodic collision chain 7 is 0{^)-shadowed 
modulo time translation by an almost collision T ^-periodic orbit of the 3 
body problem with period = T + 0{i-i). 
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Theorem 5.2. Let (s,x, € Vk be a nondegenerate modulo rotation criti- 
cal point of A^'~^ satisfying the changing direction condition US. 6\} . Then for 
sufficiently small ^ > 0, the corresponding collision chain 7 is O {11) -shadowed 
modulo rotation and time translation by an almost collision periodic modulo ro- 
tation orbit 7^ of the 3 body problem with energy E and angular momentum 
G. 

These discrete versions of Theorems I2.1H2.2I are most suitable for appHca- 
tions. In a future pubhcation we will use them in [5] to find many nontrivial 
second species solutions. 

For a dynamical systems reformulation, it is convenient to introduce Jacobi's 
discrete action functional 

n 

It is defined on 

Wif = {X = {Xi, ...,Xn)- {Xj,Xj+i) g Wy, Xn+1 = Xi}. 

Equating to the derivatives of ^^'^(s,x) with respect to s, we obtain: 

Proposition 5.4. Any nondegenerate modulo symmetry periodic collision chain 
with energy E corresponds to a nondegenerate modulo rotation critical point x 
of some J[f . 

A critical point x — {xi, . . . ,Xn) of is a n-periodic trajectory of a dis- 
crete Lagrangian system (C^) with a multivalued discrete Lagrangian = 

d 

— (if,_i {xj^i,Xj) + {xj,Xj+i)) = 0, j = l,...,n. 

Thus description of second species solutions is reduced to the dynamics of a 
discrete Lagrangian system Under a twist condition, a periodic tra- 

jectory of system corresponds to a periodic trajectory {xj,yj), yj = 

— g^L^j (xj , Xj+i), of a sequence of symplectic twist maps {xj,yj) — >■ (zj+i, j/j+i) 

with generating functions . We postpone this reformulation to a future pa- 
per, where we deal with chaotic almost collision orbits. 

The minimal degeneracy of a critical point of is at least 1 due to rotational 
symmetry J^(e*^x) = J^{x)- This implies that the discrete Lagrangian system 
(or the corresponding symplectic map) has an integral of angular momentum 
G = ixj ■ yj. One can perform Routh's reduction in this discrete Lagrangian 
system reducing it to one degree of freedom [12] , but this complicates the discrete 
Lagrangian. 

For periodic modulo rotation collision chains with fixed E, G we have: 
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Proposition 5.5. Any nondegenerate periodic modulo rotation collision chain 
with energy E corresponds to a nondegenerate modulo rotation critical point 
(x, $) of the discrete Jacobi-Routh action functional 

J^«(x, $) = J^(x) - G$, x„+i = e^*xi. 

The proofs of Theorenis l5.1l and l5.2l are modifications of the proof of Theorem 
2.1 in [5]. They are based on the Levi-Civita regularization and shadowing. The 
proof of Theorem 15.11 wih be given in the next section. The proof of Theorem 
I5.2l is similar and will be omitted. 

6 Proof of Theorem 15.11 

For /i 7^ 0, the action functional (|2.ip of system {H^) is singular when 7 ap- 
proaches A. We will formulate a variational problem for almost collision orbits 
of system {H^j) with given energy E which has no singularity at A. 

Let us fix energy E < Q. Trajectories of system (-ff^) with energy E are 
extremals of the Jacobi action functional 

J^^il) = f dsf, dsf = max{p • dq : H^{q,p) = E}. 

Away from A, the functional is a regular perturbation of the Jacobi func- 
tional for system (Hq). Regularizing J^f near A requires some preparation. 

First we describe local behavior of trajectories of system (Hq) colliding with 
A. We will use the variables (this is a version of Jacobi's variables) 

x = aiQi + 0:292, y=Pi+P2, u = q2-qi, w = aiP2 - "2^1- (6.1) 

Thus X is the center of mass of mi,m2, y is their total momentum, u is their 
relative position, and v is the scaled relative velocity. The change is symplectic: 

p ■ dq = pi • dqi + P2 ■ dq2 — y ■ dx + v ■ du. 

The inverse change is 

qi=x-a2U, qi=x + aiu, pi=aiy-v, p2 = a2y + v. 

For solutions of system (Hq), y ~ x and v = aii, where a = aia2- 

Let 7 be a trajectory of system (Hq) with energy E. We denote by {x{t), y{t),u{t), 
its representation in Jacobi's variables. If 7 has a collision ait — 0, i.e. u{0) = 0, 
a;(0) = xo, then 

Ho = \y0\y2 + \v0\y2a - \xo\~^ = E. 

We assume that collisions occurs with nonzero relative speed vq 7^ 0. Then 
there exists 5 > such that (xq, j/q) lies in a compact set 

Af = A/f -{(xo,2/o):A(a;o,yo)=i?-|yolV2+|xor^><5, \xo\>S}. (6.2) 

We fix (5 > 0. Eventually it will taken sufficiently small. Denote Bp = {u E K.^ : 
\u\ < p} and Sp = dBp. We have 
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Lemma 6.1. Take any S > and let p > be sufficiently small. Then for any 
{xo,yo) G M and any w+ G Sp there exists a trajectory 7+ : [0,r+] — >■ of 
system (Hq) with energy E such that: 

• u{t) e Bp for <t < T+ and x{0) = xq, j/(0) = yo, u{0) = 0, u{t+) = u+. 

• 7+ smoothly depends on {xo,yo,u+) G M x Sp and 

T+ = T+{xo,yo,u+) = p^/a/2X{xo,yo) + 0{p'^), 
x(r+) = £,+ {xo, yo, w+) = xq + p^/a/2X{xo,yo)yo + 0{p'^). (6.3) 

• The Maupertuis action of 7+ has the form 

-^^(7+)=/ P- dq = a+{xo,yo,u+) (6.4) 
"'7+ 

= p^j2a/X{xo,yo){E + \xo\-^) + 0{p^). 

Remark 6.1. On Sp we use the polar coordinate 6, where u = pe'^. Thus 
OIp^) m,eans a function of XQ,yo,0 whose norm is bounded by cp^ with c 

independent of p. 

The proof is obtained by a simple shooting argument, because Hq has no 
singularity at A: 

x{t) =xo + tyo + 0{t^), u{t) = tvo/a + 0{t^). 

It remains to solve the equation w(t+) = u+ for t+ and vq, where |fo| = 
^2aA(a;o,yo)- □ 

Similarly, wc have a trajectory 7_ : [r_, 0] U"^ of system (Ho) with energy 
E such that a;(0) = xq, y(0) = yo, m(0) — 0, u(r_) — u-. Then 

r_ = T_(a;o,yo,M_), x{t-) = (,-{xo,yo,U-), J^(7_) = a_(xo, t/o, "-)• 

If (xcj/o) G M, then xo belongs to 

D = Df = {x:5<\x\<{5- E)-^}. 

Let be the boundary of the tubular neighborhood Np of C A: 

T.p = {q: X € D, it G Sp], Np = {q : x & D, ue Bp}. 

Fix arbitrary large C > and let 

Kp = {iqo,q+) G X : |(7o - g+l < Cp}. 

Lemma 6.2. If p > is sufficiently small, then for any (q^.q^) G Kp, there 
exists a trajectory 7 : [0, T-^.] — )■ Np of system (Hq) with energy E joining qo with 
q+. Moreover 'y smoothly depends on {qo,q+) and its Maupertuis action has the 
form 

J^ii) = dE{qQ,q+) 
= y/2a-^{E + |xo|-i)(|x+ - xo|2 + a2p2) + ^(p^). (6.5) 
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Here dE{qoi<l+) is the distance in the Jacobi metric dsE- 

Proof. The condition 7(t-|_) = = (x+j gives u{t^) = x(r+) = 
Then (|6.3p makes it possible to determine 



2/0 = V+i<lo,q+) = \ I ] - xo) + 0{p). (6.6) 

y \x+ - xor + a P 

□ 

Next we connect points q-,q+ G Sp by a reflection trajectory of energy E. 
Let Pp = {(<z_,q+) e : |g_ - 9+] < Cp}. 

Proposition 6.1. Let p > be sufficiently small. Then for any {q^,q^) G Pp: 

• There exist r_ < < t_|_ and a broken trajectory 7 : [r_,T-^] — > A'p 
wii/i energy E such that 7(0) = go = (2^0,^0) € -D^, 7|[o,r+]7 7|[t_,o] '^'^g 
trajectories of system (Hq), 7(t±) = q± and there is no jump of total 
momentum at collision: y{+0) = y(— 0) = yo. 

• 7 smoothly depends on {q-,q+) G Pp. 

• The Maupertuis action of 7 has the form 

JE{n)= dsE ^ gE{q-,q+) ^ dEiqQ,q-) + dEiqo:q+) (6.7) 



= V2{E + \xo\~^){\x+ - x^\^ + Aa^p^) + 0{p^). (6.8) 

xo = aq-,q+) = ix+ + x-)/2 + o{p^), (6.9) 



l2{E+\xo\-^) ^^n(^ 



Proof. We find from the equation 

— {dE{qo,q-)+dE{qQ,q+)) = ^ ?7+(qo, g+) = ??- (go, <?-), 
where 77-1- is defined in (|6.6p . Differentiating (|6.5p , we see that the Hessian matrix 



8{E+\xo\-^) [x+~x^)®{x+-x^) ^ ^^^^ 



\j a{\x+ — x^\'^ + Aa^ p^) \ |a;+ — p + 4a2p2 

is nondegenerate. By the implicit function theorem, the solution xo = $(<?-, q+) 
is smooth. □ 

A similar result holds for the perturbed system {Hp), but it is no longer easy 
to prove. Fix an arbitrary small constant 5 > Q. 
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Lemma 6.3. Let p > be sufficiently small. There exists fiQ > such that for 
all ji e (0,/io], any {xo,yo) £ M and any u± £ Sp such that + u_| > Sp: 

• There exist t_ < < i+ and a trajectory 7 : [t_, — >■ Np of system [Hp) 
with energy E such that u{t±) = u±, x{0) = xq, y(0) = j/o- 

• 7 smoothly depends on {xg, yQ,u^,u^, p,) G M x SpX (0, po] and converges 
to a concatenation of trajectories "f± in Lemma \6.1\ as p 0. 

• The Maupertuis action of 7 has the form 

J^{l)= / P - dq ^ fpixQ,yQ,u.,u+) (6.10) 
= a+{xo,yo,u+) + a^{xo,yo,u^) + pa{xo,yo,U-,u+, p), (6.11) 
where a is hounded on AI x Sp x (0, pq]. 

t± = T±{xo,yo,u±) + pf±{xo,yo,u+,u^,p), 
x{t±) = x^{xQ,yQ,u_,u+) 
= £,±{xo,yo,u±) + px±{xo,yo,u-,u+, p). (6.12) 

where t± and x± are uniformly bounded on M x x (0,/uo]. 

• // \u+ — u_| > Sp, then 

pa< min dH(t),A)<pb, < a < b. (6.13) 

tG[t_,t+] 

The proof of Lemma 16.31 is given in section [T] It is based on Levi-Civita 
regularization and a generalization of Shilnikov's Lemma [20], see also [H], to 
normally hyperbolic critical manifolds of a Hamiltonian system. 

Next we deduce a local connection theorem. Fix arbitrary small 6 > 0, 
arbitrary large C > and let 

Qp = ■■ \q- - q+\ < Cp, \u+ + u_| > 5p}. 

Theorem 6.1. Let p > be sufficiently small. There exists /ig > such that 
for all {q^,q+,p) eQpX (0,/io].- 

• There exist <_ < < and a trajectory 7 : [t_, i+J — > Np of system (Hp) 
with energy E such that 7(i±) = q± and the minimum of d{j{t), A) is 
attained at t — 0. 

• 7 smoothly depends on (g-, q+, p) € Qp x (0, po] and converges to a reflec- 
tion trajectory in Provosition \6. 1\ as p ^ 0. 
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• The Maupertuis action of 7 has the form 

•^/f (7) = / P ■ dq ^ gjfiq-,q+) ^ gE{q-,q+) + tJ-9iq-,q+,t^), 

where g is hounded on Qp x (0, ^0] . 

• // |u+ — U-l > Sp, then i6.13\) holds. 

Thus the action 0/^(7) — gp{q-,q+) has a hmit gE{q-,q+) as fi which 
is the action of the reflection orbit in Proposition 16.11 The condition that the 
distance to A is attained at i = is needed only to exclude time translations, 
so that t± are uniquely defined. 

Proof. We need to find {xo,yo) such that xj^(xo, yoi ^-i "+) Since the 

implicit function theorem worked in the proof of Proposition 16.11 by (j6.12l) , for 
small ^ > it will work also here. □ 

Proof of Theorem I5.il Let 7 be a nondegenerate n-collision chain with energy 
E. Let t = (^1, . . . , t„) be collision times, x — (xi, . . . , a;„), j{tj) = {xj,Xj), the 
corresponding collision points and Uj the collision total momenta. Take S > 
so small that the collision points and the collision speeds = v{tj ± 0) satisfy 

\xj\>S, \v+\^\vJ\>SV2^. 

Then {xj, yj) £ M and Xj £ D. 

Take small p > and let = tj ± be the closest to tj times when q^ = 
litf) £ 5]p. Since 7 satisfies the changing direction condition, {q~,q^) £ Qp 
if C > is taken sufficiently large and p > sufficiently small. Moreover for 

£ Ep close to q^, we have £ Qp. Thus by Theorem 16.11 for small 

fj, > the points can be joined in Np by a trajectory of system (iJ^) with 
energy E and the Maupertuis action J^(7^) = 

Since ^\[tj,tj+i] is nondegenerate and, by no early collisions condition, does 
not come near A, for ^+ close to and close to q~_^^ and small p > 0, the 
points ^j' and can be joined by a trajectory cr^ of system (Hp) with energy 
E and the Maupertuis action Jp[<Jp) = ^^("^j^, Cj+i)- This trajectory smoothly 
depends on /i also for /i = 0, and ^o" (Cj^' Cj+i) the Maupertuis action of a 
connecting trajectory of system {Hq). 

Combine the trajectories 7^^, ct^ in a broken trajectory 7^ with energy E 
and Maupertuis action 

n 
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The function has a hmit /o as /i — > and 
where 

n 

and = ^(^~,^^) is defined by (|6.9p . The remainder f{£,,n) is bounded on 
y X (0,/io], where the neighborhood Y C S^" of 

is independent of ^. 

Looking for critical points with respect to with fixed £,j we obtain /o = 
J]f (^1, . . . , for some k e Z^", and Jjf has a nondegenerate modulo rotation 
critical point x. Thus /o has a nondegenerate critical point q S S^". Then for 
small /i > the function /p(^) has a nondegenerate modulo rotation critical 
point close to q. The corresponding broken trajectory 7^ has no break of 
velocity at intersection points with Ep and hence 7^ is a periodic trajectory 
of system (iJ^) with energy £'. □ 



7 Levi-Civita regularization 



In this section we prove Lemma [^751 In the Jacobi variables (|6.ip . the Hamilto- 
nian 77^ takes the form 

^ _ (1 I bp «i «2 f^a 

^ 2 2a \a2U — x\ \aiu + x\ \u\ 

Let us perform the Levi-Civita regularization on the fixed energy level i/^ = 
E. We identify u, f G = C with complex numbers and make a change of 
variables 

u = f, v = ii/2l 

Since 

V ■ du = Re (v du) = Re (77 0?^) = -q ■ d^, 
the change is symplectic: 

p ■ dq — y ■ dx + T] ■ d^. (7-1) 
Finally, we obtain a transformation 

9i=2^-a2^^, qi^x + aiS,'^, pi ^ aiy - r]/2^, p2 = "2^ + ??/2C- 
The Levi-Civita map 

g : R2 X M2 X [/ X M2 ^ (R4 \ A) X R\ g{x,y,i,rj) = (91,92,^1,^2), 
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is a symplectic double covering undefined at ^ = which corresponds to the 
collision set A. 
Let 

H'',{x,y,^,Jl) = \£,\^{H^og-E)+^ia 

Let = {i/^ = £;} and L^ = {H^ = fia}. Since ^(r^) = S^, the map 
g takes orbits of system (H^) on to orbits of system (-ff^) on S^. The 
time parametrization is changed: the new time is given by dr = In the 

following we will continue to denote the new time by t. 

The singularity at A disappeared after regularization. The regularized Hamil- 
tonian is smooth on 

which means excluding collisions of mi and mi with m^. The parameter /xa 
may be regarded as new energy. The rotation group and the integral of angular 
momentum are now 

{x, y, 77) ^ {e^'x, e'y, e'^/^^, e^^/^yy), G = • j/ + ^^ 77/2. 

The Hamiltonian T-L^ has a critical manifold ^ = 77 = which is contained 
in the level set of 71^ . We have 

2/, ^, = ^ - iepA(x, 2/) + Om^). 
8a 

Collisions of ttii , 7712 with nonzero relative velocity correspond to the solutions 
asymptotic to 

X = M X {(0,0)}, 

where M is as in (j6.2p . This is is a compact normally hyperbolic symplectic 
critical manifold for Ti,^ . We obtain 

Theorem 7.1. Collision orbits of system (Hq) with energy E correspond to 
orbits of system {Hq) doubly asymptotic to A4. Orbits of system {H^) with 
energy E passing 0{ii)-close to the singular set A correspond to orbits of system 
(H^) on the level F^ passing 0{^JJi)-close to A4. 

Next we translate Lemma [6.31 to the new variables. 

Let r > and let Afr = M x he a tubular neighborhood of in T'. 
By the stable and unstable manifold theorems for normally hyperbolic invariant 
manifolds [13], if r > is small enough, for any {xo,yo) G M and ^_ G Sr there 
exists a unique solution C- : [0,-|-cx)) -> A/'r, C-{t) — {x{t),y{t),^{t),r]{t)), of 
system {Hq) such that ^(0) — ^_ and (^{00) — {xQ,yQ, 0, 0) & M. We denote its 
action by 

J{C-)^J y ■ dx + 7] ■ d£_ ^ J-{xo,yo,£_^). 
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Since the stable and unstable manifolds are smooth, J_ is a smooth function 
on M X Sr- Similarly we define the function J+ on M x Sr as the action of a 
solution ^+ asymptotic to as i — > — cx). 

We have an analog of Shilnikov's Lemma j20j . Fix small e > and denote 

Qr = {(xo,yo,^-,C+) e M X S?. : ^- ■ ^+ > eV}. 

Theorem 7.2. There exists r > and /ig > suc/i that for any {xq, t/Oi C-i C+j /^) ^ 
Qr X (0, /lo] •■ 

• There exists T > and a solution 

at) = {x{t), y{t),at), V{t)) eK, t€ [~T, T], 
of system (H^) on such that 

x{0)=xo, 2/(0) -yo, a-T) = ^-, e(T)=e+. (7.2) 

• C smoothly depends on (xq, j/o, m) G Qr x (0,/io]- 

• The Maupertuis action is a smooth function on Qr x (0,/io] o,nd has the 
form 

■HO ^ J y-dx+iydi = J^{xo,yQ,i-)+J+{xQ,yQ,i+)+^iJ{xo,yQ,i-,i+,^i), 

(7.3) 

where J is hounded on Qr x (0, /ig] • 

A result very similar to Theorem 17.21 was proved in jB]. A complete proof of 
Theorem 17.21 will be published in [5]. 

Proof of Lemma WT^ We set p — r^ and u = For given u± G Sp take ^± e S'r 
such that ^+ • ^_ > 0. If ^+ • > 0, then u+ 7^ — u_. Moreover for given S > 
there exists £ > such that \u- + u+\ > Sp implies ^_ • ^+ > e^r^. If ( is 
a trajectory in Theorem 17. 2[ then the corresponding trajectory (/(^) of system 
(iJp) satisfies the conditions of Lemma 16.31 In particular, by (|7.1I) . 

•/±(a;o,2/o,C±) = a±(a;o,2/o,'"±)- 

□ 
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